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Rotational invariance is a well-established feature of low-energy physics. Violations of this sym- 
metry must be extremely small today, but could have been larger in earlier epochs. In this paper 
we examine the consequences of a small breaking of rotational invariance during the inflationary 
era when the primordial density fluctuations were generated. Assuming that a fixed-norm vector 
picked out a preferred direction during the inflationary era, we explore the imprint it would leave 
on the cosmic microwave background anisotropy, and provide explicit formulas for the expected 
amplitudes {ai m a*i m i) of the spherical-harmonic coefficients. We suggest that it is natural to expect 
that the imprint on the primordial power spectrum of a preferred spatial direction is approximately 
scale- invariant, and examine a simple model in which this is true. 



I. INTRODUCTION 

Inflationary cosmology, originally proposed as a so- 
lution to the horizon, flatness, and monopole problems 
[H 0> provides a very successful mechanism for gen- 
erating primordial density perturbations. During in- 
flation, quantum vacuum fluctuations in a light scalar 
field are redshiftcd far outside the Hubble radius, im- 
printing an approximately scale-invariant spectrum of 
classical density perturbations || H|. Models that re- 
alize this scenario have been widely discussed @, H, 0]- 
The resulting perturbations give rise to galaxy forma- 
tion and temperature anisotropies in the cosmic mi- 
crowave background, in excellent agreement with obser- 
vation @, S M, EE S3, 0, 0, 0, um ■ 

If density perturbations do arise from inflation, they 
provide a unique window on physics at otherwise inac- 
cessible energy scales. In a typical inflationary model 
(although certainly not in all of them), the energy scale 
E = V 1 / 4 is related to the amplitude of density fluctua- 
tions 5 and the reduced Planck mass Mp via E ~ VSMp. 
Since we observe 5 ~ 10~ 5 , it is very plausible that infla- 
tion occurs near the scale of grand unification, and not 
too far from scales where quantum gravity is relevant. 
Since direct experimental probes provide very few con- 
straints on physics at such energies, it makes sense to 
be open-minded about what might happen during the 
inflationary era. 

In this paper we ask what happens when a cherished 
property of low-energy physics - rotational invariance - 
is violated during inflation. Rotational invariance is of 
course a subset of Lorentz invariance, and theoretical 
models of Lorentz violation in the current universe (and 
experimental constraints thereon) have been extensively 
studied in recent years [13, 0, 0, 0] ■ Here we are 
specifically concerned with the possibility that rotational 
invariance may have been broken during inflation by an 
effect that has subsequently disappeared, and study the 
effects of such breaking on CMB anisotropies. It is possi- 
ble that such an effect has already been detected, in the 
form of the "Axis of Evil," an apparent alignment of the 
CMB multipoles on very large scales [H 0, 0, 0, 0, 



0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, m . 
Although its statistical significance is hard to quantify, a 
variety of models have been pu t forward to ex plai n this 
phenomenon [H 0, 0, 0, 0, 0, 0, 0, 0,13 . Our 
aim is not to construct a model contrived to explain the 
currently observed large-scale anomalies, but rather to 
make robust predictions for the observable consequences 
of a preferred direction during inflation, allowing obser- 
vations to put constraints on its magnitude. 

The power spectrum P(k) for the primordial density 
perturbations 5(k) is defined by 

(<5(k)r(q)>=P(fc) ( 5 3 (k-q). (1) 

The Dirac delta function in Eq. ([1]) implies that modes 
with different wavenumbers are uncoupled, and is a con- 
sequence of translational invariance during the inflation- 
ary era. On the other hand, the fact that the power 
spectrum P(k) only depends on the magnitude of the 
vector k is a consequence of rotational invariance. Sup- 
pose that during the inflationary era rotational invariance 
is broken by the presence of a small vector that points 
in the direction of a unit vector n. Assuming a parity 
k — > — k symmetry, the leading effect of the violation of 
rotational invariance changes the most general form of 
the power spectrum from P(k) to P'(k), where 

P>(k)=P(k)(l + g(k)(k-nf). (2) 

Here k is the unit vector along the direction of k and 
we are neglecting higher powers of k • n since they will 
be suppressed by more powers of the magnitude of the 
small vector that breaks rotational invariance. (Effects 
of a timelike vector on inflationary perturbations have 
also been studied [Hj].) 

Towards the end of the inflationary era, the physical 
wavelengths that correspond to scales of astrophysical 
interest are large compared with the inverse Hubble con- 
stant during inflation or any of the dimensionful particle- 
physics quantities that might be relevant during inflation. 
The same naturalness arguments that lead to the scale- 
invariant Harrison- Zeldovich spectrum (i.e., primordial 
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P(k) oc 1/fc 3 ) imply that g(k) in Eq. (0) should be in- 
dependent of k. Assuming that g(k) is a fc-independent 
constant g» over the scales of astrophysical interest, we 
arrive at 



P'(k) =P(fc)(l+. 9 ,(k-n) : 



(3) 



This is the form of the primordial power spectrum that 
takes into account the leading effects of the violation of 
rotational invariance by a small vector in the inflationary 
era that points in the direction n. In the next section we 
discuss the implications of the power spectrum in Eqs. @ 
and ([3]) for the anisotropy of the microwave background 
radiation. The breaking of rotational invariance gives rise 
to correlations between multipole moments that would 
normally vanish and also alters the predictions for the 
usual multipole moment correlations. In section three 
we discuss a simple model that realizes the form of the 
primordial power spectrum in Eq. ([3]). Concluding re- 
marks are given in section four. 



II. MICROWAVE BACKGROUND 

We are interested in a quantitative understanding of 
how the substitution, P(k) — > -P'(k), changes the pre- 
diction for the microwave background anisotropy AT/T. 
The multipole moments are defined by 



air. 



J dn.(ir(e))*Y(e) 



(4) 



The anisotropy of the microwave background tempera- 
ture T along the direction of the unit vector e is related 
to the primordial fluctuations by 

^(e) = J dk£ (H!±i) {-i) l P l {^.e)5{\ S )@ l {kl 

where Pi is the Legendre polynomial of order I and 0/ (fe) 
is a function of the magnitude of the wave- vector k that 
includes, for example, the effects of the transfer function. 
It can only depend on the magnitude of the wave- vector 
since the dynamics after the inflationary era is assumed 
to be rotationally invariant. 

We would like to compute the expectation values 
( a im,a*/ m i) to first order in the small quantity g(k) that 
characterizes the primordial violation of rotational invari- 
ance. We write 



{aimalm') = (oj m o*/ m />o + A(7m; I'm'), 



where the subscript denotes the usual rotationally in- 
variant piece, 



(oj m o*/ m /)o = Sii'S- 



,,..„> i >\kk-n<k^i<kr. (7) 

o 



It is useful to introduce the "spherical" components of 
the unit vector n that defines the preferred direction for 
rotational non-invariance, 



n + 



V2 



V2 



, n = n z . 



. . (8) 

In terms of these components the unit norm condition 
becomes — 2n + n_ = 1. Note that we do not assume 
that the preferred direction n coincides with the z axis 
of the coordinate system used to parameterize the mi- 
crowave sky (i.e., that n + = n_ = 0). Expressions anal- 
ogous to ours have been derived by Gumrukcuoglu et al. 
|52| under the assumption that these two directions are 
coincident; see also [37| . 

Using the identity 



4-7T 

m— — I 



(9) 



it is straightforward to express the sought-after pertur- 
bation as 



A(!m;!W) = (-f'e, m;W / dM 2 P(fc)<?(fc)eKfc)@;'(*0> 

Jo 

(10) 

where 



Ait 

Zlm;l>m' = y / d^F^k))*^ 1 (k) (III 



n+y^k) + n-Yf x (k) + n Y?(k) 



(we use the Condon-Shortley phase convention for the 
spherical harmonics; see [5J]). 

The integral in pop encodes information about the 
power spectrum and the transfer function, as well as the 
scale-dependence of the preferred-direction effect, while 
the constants £,i m ;Vm' are purely geometric. The integra- 
tion over solid angles is straightforward to perform. It is 
convenient to decompose the £i m - t i' m ' into coefficients of 
the quadratic quantities riiTij, via 



(6) ilm-l'm' 



+ »-€,-> m / + 2n+n_£+- rm , (12) 



n 2 f ++ 

"■+ ( >lm;l'm' 



These coefficients are then given by the following expres- 
sions: 
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^lm\l'm' 



f ++ 

Hmd' 



= -8 m ' 



slmd'm' 



m',m+2 



^/(P - (m + l) 2 )(Z + m + 2)(Z-m) 1 / (Z + m + 1)(Z + m + 2) (/ + m + 3)(Z + m + 4) 

' V (2I + 3)(2I-1) 2 l '' l+2 V (2Z + 1)(2Z + 3) 2 (2Z + 5) 



1 (I - m)(l - m - - m - 2)(l - m - 3) 

2 l '' l ~ 2 V (2Z + 1)(2Z - 1) 2 (2Z - 3) 



I' rn':lm ' 



or (-1 + Z + Z 2 + to 2 ) . ((Z + l) 2 -m 2 )((/ + 2) 2 -m 2 ) 

-^i',* - T^Ti ,n/o, , on 1" °l',l+2\ 



(2Z — 1)(2Z + 3) 



(2/ + 1)(2Z + 3) 2 (2/ + 5) 



f (l 2 -m 2 )((Z-l) 2 -ro 2 ) 
(2Z — 3)(2Z — 1) 2 (2Z + 1) 



. {2m + l)J(l + m + l)[l-m) . / ((Z + l) 2 - m 2 )(Z + m + 2){l + m + 3) 

Oi'.i 7771 , on ^ <V,!+2W — 



(2Z — 1)(2Z + 3) 



(2Z + 1)(2Z + 3) 2 (2Z + 5) 



8l'l-2\ 



1 (I 2 - m 2 )(l - m- 1)(Z - m- 2) 
(2/ — 3)(2/ — 1)2(2/ + 1) 



~ Si'm' 



V mJ'.Vm ' 



Hmd'm' ~ u m,m' 



. (2l 2 + 21 - 2m 2 - 1) . /((/ + l) 2 -m 2 )((/ + 2) 2 -m 2 ) 

d U' 77^1 ,n/^, , on r-(ty,J+24 



(2Z-l)(2Z + 3) 



(2Z + l)(2Z + 3) 2 (2Z + 5) 



+8l'l-2\ 



(I 2 -m 2 )((Z -l) 2 -m 2 ) 
(2/ - 3)(2/ - 1) 2 (2/ + 1)) 



(13) 



The formulas (j 1 21 1 3[) are explicit expressions for the ge- 
ometrical part of the perturbation (|10|) . As we mentioned 
in the introduction, it is natural to imagine that the vio- 
lation of rotational invariance is approximately scale in- 
variant, which implies that it is a good approximation to 
set g{k) = g*, a constant. If we define polar coordinates 
9*,(f>* for the preferred direction, 



sinZLcoso 



sinZLsino 



cos6>* 



(14) 



these expressions can be compared directly with obser- 
vations to constrain the three parameters (g* , 6>* , 0* ) . 

When g(k) — g*, a simplification occurs for / = /' and 
m = mf, as the dependence on the power spectrum for 
the terms that violate rotational invariance A(/m; Im) 
is the same as the rotationally-invariant part (a; m a ; * m )o- 
We can then find a simple expression for their ratio, 



sin 



A(/m;/m) g* 
{ai m a*i m )o 2 

(3cos 2 Z 



(15) 



2Z 2 + 2Z - 2m 2 - 1 
(2/ — 1)(2/ + 3) 



For large multipoles, Z 3> 1, and for the magnitude of m 



of the order of Z, this expression simplifies to 



A(lm;lm) g* 
(ai m a* m ) 4 



1 + cos 



(3cos 



. in 



(16) 



III. INFLATION MODEL WITH A PREFERRED 
DIRECTION 

It is interesting to see how the rotationally non- 
invariant power spectrum in Eq. ([3]) can arise in an ex- 
plicit model of anisotropic inflation. We will assume that, 
during most of the inflationary era, rotational invariance 
is broken by a spacelike four-vector u M with invariant 
length 



g llv u iL u v = m 2 . 



(17) 



We will consider the effect of the energy-momentum ten- 
sor associated with this vector on the expansion of the 
universe during inflation, ignoring direct couplings of 
w M to other fields. Gravitational effects of dynamical 
Lorentz-violating vector fields have been considered pre- 
viously in the literature [H, H3, IH IH IM H3 • 

We assume that the four-vector is non-zero only 
during the time interval < t < /*, where i* is the end 
of inflation, so that the dynamics is rotationally invariant 
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during reheating and thereafter. During the time interval 
< t < t*, the dynamics of interest is governed by the 
action 



where 



and 



(18) 



(19) 



-0x Wu a V M iv - 2 (V^) 2 (20) 
-#,V"u a V ff u M + - m 2 ) . 



Here A is a Lagrange multiplier that enforces the con- 
straint (fT7|) . Quantum fluctuations in the massless scalar 
field x are assumed to dominate the density perturba- 
tions via the DGZK mechanism [55j |. In that case we 
need simply calculate the fluctuations in x, without wor- 
rying about the behavior of the inflationary potential. 

We approximate the inflaton energy density as a con- 
stant, modeling the effects of the inflaton field by a vac- 
uum energy p\ in Eq. (|18p . The inflationary spacetime 
is taken of the form 



ds 2 



-dt 2 + a(t) 2 dx 2 L + b(t) 2 dz 2 



(21) 



since we have chosen the four- vector to be aligned along 
the z-axis direction, 



u° = , u x = 



,u z 



m (22) 



The energy-momentum tensor for u M derived from (|20p 
is [11] 

T$ = 2/3i(V p u'V I ,u„-Vu / *V < ,u„) 

-2[V p (w (At J' , , ) ) + V P (/J (H ) - V p (u (M J,) p )] 
+2m~ 2 u,jV pJf'u^Uv + g^C u , (23) 

where J^o- is the current tensor, 

= -A W - ft <S£ V p u>> - foVru". (24) 

Given (|22p and ((23J), the nonvanishing components of 
the stress tensor are 



= Tyy = P im 2 a 2 



m iU 2 -2bb-4^ 



(25) 



Note that the components of the energy momentum ten- 
sor in our chosen background are independent of 02 and 

03- 



Solving Einstein's equation during the time interval 
< t < t*, with initial conditions a(0) = 1 and 6(0) = 1, 
gives 



o(t) = 



&(*) = e Wbt , 



(26) 



where 



H a = - = ff 6 (l + 16ttG/3ito 2 ), 
a 



As = 



8irGp\ 



b V (1 + 87tG/3ito 2 )(3 + 327rG/3im 2 ) 



.(27) 



According to the cosmic no-hair theorem, initially ex- 
panding homogeneous cosmological models in the pres- 
ence of a positive cosmological constant will rapidly ap- 
proach a de Sitter solution, if the other matter fields obey 
the dominant and strong energy conditions [62j . Our spe- 
cific model violates these conditions. Nevertheless, for 
03 = —0i and 02 = the kinetic term for fluctuations 
about our background has the form of a field strength 
tensor squared and so is ghost free. We therefore ex- 
pect the configuration to be stable with respect to small 
fluctuations. 

It will turn out to be convenient to refer to a fictitious 
isotropic metric, 



ds 2 



-dt 2 + a(t) 2 [dx 2 + dy 2 + dz 2 }, 



(28) 



H = ^(2H a + H b ). 



in which the scale factor expands exponentially 

o(t) = e m (29) 
with an "average" Hubble parameter, 

(30) 

Deviations from isotropy can be parameterized by 

2 / H b - H a \ 
eff= 3l v T _ J' (31) 

where the 2/3 will become useful later. We work in the 
limit iV*|efj| << 1, where N* = Ht* is the number of 
e-foldings during the time when the four-vector is 
non-zero. This assures that the violation of rotational 
invariance due to the anisotropic expansion is always a 
small perturbation. 

We need to compute the correlation function 
(x(x, t)x{y, t)). Treating ch as a small perturbation, we 
find that to first order in this quantity we obtain (see for 
example ref. [63l |) 



(x(x,t)x(y,t)> - (xi(x,t)xi(y.*)> 



(32) 



-* / dt'([iJ / (t'),x/(x,t)x/(y,i)])- 
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Here the interaction-picture Hamiltonian Hi(t) is given 
by 



Hi{t) = J d?x\ 



ait) 



\dx 3 



(33) 



The interaction-picture (i.e. free) field obeys the 
rotationally-invariant equation of motion, 



d 2 Xi 
dt 2 



3H 



dt d(t) 2 dx 2 



= 0. 



(34) 



We can write the two-point correlation function |33|) 
in terms of Fourier transforms as 



(x(x,t)x(y,t)) 



d 3 fc 



-tk-(x-y) 



(2tt) 3 

(k • n) 2 AP(fc) 



P(k) 



(35) 



Converting to the conformal time of the isotropic metric, 

1 



-Ht 



H 



(36) 



and expanding in e#, we find that P(k) ~ Ixi ^ 7 ")! 2 ) an d 



AP(k)~3ik 2 e H f dr'f--^ 



(37) 



xlog(-tfT') { X f{r')xT (T)*f (xr(rrxt"(r)) 



(0) 



,(0), 



where 



Xi 0) (r) 



— ikr 



2k 



T -~k 



(38) 



We assume that the modes k of astrophysical interest 
have wavelenghts much smaller than the Hubble radius 
at the beginning of inflation, which in our normalization 
implies k >> H. They cross the horizon around sixty 
e-foldings before the end of inflation (which we take to 
occur at about £*). Taking |fcr| << 1, we find that 



AP(k) ~ ^e H ^log(k/H), 



(39) 



where we have neglected contributions not enhanced by 
the large logarithm. 

There is another way to derive Eq. ([39]) . For modes 
with wavenumbers along the z direction or perpendicular 
to this direction, the Fourier transform of the two point 
function (x(x, t)x(y, t)) can be found exactly without re- 
sorting to perturbation theory. For example, modes Xk 
with k = kz (wavevectors parallel to the preferred direc- 
tion) obey the differential equation 



d 2 Xk 
dt 2 



- dxk k 2 
3#-^ + TTT^Xk = 



dt 



b{tf 



(40) 



The canonical commutation relations imply that Xk sat- 
isfies the normalization condition, 

Xk(r) - — t~ — XkiV = -i(Ht) . 



\ dr 



V dr 



(41) 

We find that the properly normalized solution to Eq. (l4"01) 
is 



Xk{r) 



.3/2 



(k/H)' eH {kT) 
l + e H 



1+6H 



(2) 

where H„ is a Hankel function, and 



2 + 2e H 



, (42) 



(43) 



The contribution to the Fourier transform of the two 
point x correlation for a mode along the z direction is 
|X/c(t)| 2 - For small en and |fcr| and large k/H, this be- 
comes 



\Xk(rW 



2k 3 



1 + 3e H \og(k/H)) . 



(44) 



Here we have neglected terms linear in en that are not 
enhanced by the large logarithm. Combining this result 
with a similar analysis for modes perpendicular to the z 
direction reproduces the result in Eq. (f3"9"|) . 

Finally we note that the density perturbation power 
spectrum is defined by a Fourier transform with respect 
to coordinates where physical laws have manifest rota- 
tional invariance. However at time t = the coordi- 
nates in Eq. (|2ip do not exhibit manifest rotational in- 
variance, due to the difference between a(i») and b(t„). 
Rescaling coordinates, z — > z(a(t*)/6(t*)) and xi — > 
xj_(a(i*)/a(i*)), we find that the function g(k) charac- 
terizing the rotationally non-invariant part of the power 
spectrum for the primordial density perturbations is 



g(k) = -e H (log(k/S)-N* 



(45) 



where the term proportional to comes from the rescal- 
ing of coordinates and q(t*) = fc/a(t*) is the physical 
wavelength of the mode of interest at the end of infla- 
tion. 

The logarithm in (|45[) is actually nearly constant over 
values of q(t*) of astrophysical interest. The range of 
q{t*) probed by CMB measurements is about a factor 
of 10 3 , so log(q(t*)/S) changes by roughly 7. But the 
modes of cosmological interest cross the deSitter hori- 
zon around 60 e-foldings before the end of inflation. So 
| \og{q{t*) / H)\ is approximately 60. Hence, in this model 
g(k) varies by about 10% over the range of modes of 
cosmological interest and our general expectation that 
setting g(k) = g* is a reasonable approximation has been 
confirmed. 
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For simplicity in this analysis we neglected terms that 
directly couple w M to \- For example we could have added 
the term u^u" d^x^uXI 'M 2 to the Lagrange density. It is 
easy to see that this gives an additional scale invariant 
contribution, 3m 2 /M 2 , to g(k). 

IV. CONCLUDING REMARKS 

We have investigated the possibility that rotational in- 
variance may have been explicitly broken during inflation 
by an effect that has disappeared in the later universe. 
The observed CMB temperature anisotropics provide a 
direct window onto the physics of the inflationary era, 
and therefore offer a unique opportunity for constraining 
(and discovering) new phenomena at high scales. Our 
aim has been to investigate the generic predictions we 
expect from the presence of a preferred direction during 
inflation. 

If rotational invariance is violated during inflation, it 
is natural for the effects of such a violation to show up 



in a scale-invariant way, just as the amplitude of the per- 
turbations themselves are approximately scale- invariant. 
Under that assumption, we derive a powerful set of pre- 
dictions for the expectation values (ai m a*, m ,) that de- 
pend on only three parameters: a single amplitude g*, 
and a direction on the sky defined by a unit vector n. In- 
vestigation of a simple model confirms the approximate 
scale-independence of this effect. The resulting expres- 
sions (|10I12I13|) can be directly compared with observa- 
tions to probe the existence of small Lorentz-violating 
effects in the very early universe. 
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